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The light scattering spectra of molecular liquids are derived within a generalized hydrodynamics. 
The wave vector and scattering angle dependences are given in the most general case and the change 
of the spectral features from liquid to solidlike is discussed without phenomenological model assump- 
tions for (general) dielectric systems without long-ranged order. Exact microscopic expressions are 
derived for the frequency-dependent transport kernels, generalized thermodynamic derivatives and 
the background spectra. 



PACS numbers: 78.35.+C, 64.70.Pf 
I. INTRODUCTION 

Light scattering is a powerful tool to study the dynam- 
ics of dense (transparent) materials M. The fluctuations 
of the dielectric tensor for wave vector transfer q are 
measured, where the corresponding wavelength can be 
considered large compared to molecular length scales. In 
this case the theoretical description of the light spectra 
simplifies as it suffices to determine the lowest orders in 
wave vector q only. The wavevector and scattering-angle 
dependence of the scattering cross sections thus can be 
determined. 

Focusing on low frequencies, the hydrodynamic ap- 
proach can be used to calculate the spectral shapes. 
For the polarized light scattering spectra of gases and 
liquids, the assumption that dielectric fluctuations are 
dominantly caused by density fluctuations leads to the 
well known Rayleigh-Brillouin spectra. The correspond- 
ing hydrodynamic depolarized spectra were first obtained 
within simplified models by Andersen and Pecora [^) and 
Keyes and Kivelson j3|. In the latter case, upon the 
assumption that fluctuations in the orientations of the 
molecules cause the dielectric variations. As observed 
experimentally in molecular liquids, a negative central 
line appears called "Rytov-dip" |4|, |) . 

The spectra, in particular the depolarized ones, change 
qualitatively if the relaxation times of the structural dy- 
namics increase upon cooling the liquids. For the po- 
larized spectra Mountain introduced a frequency depen- 
dent longitudinal viscosity in order to model the addi- 
tionally appearing central line |(|. In the depolarized 
light scattering spectra transverse sound peaks become 
visible as expected from hydrodynamic calculations for 
solids 0. Within phenomenological models including 
non-hydrodynamic variables, the changes of the spec- 
tra from liquidlike at high to solidlike at low tempera- 
tures could be explained p|-p^[. However, assumptions 
about the included slow variables and about phenomeno- 
logical kinetic equations coupling their time dependences 
were required. The introduction of memory functions, in 



most detail in the recent work of Dreyfus et al. ^3) , 
has relaxed the requirement to identify and include all 
slow variables but still uses phenomenological equations 
to couple rotational and translational degrees of freedom. 

Here, we clarify that generalized hydrodynamics and 
symmetry considerations suffice to explain the light scat- 
tering spectra qualitatively as they change from liquid- 
to solidlike. The results will be derived without any as- 
sumptions about non-hydrodynamic variables (and their 
couplings) and will also not depend on specific light scat- 
tering mechanisms nor molecular parameters like shape, 
dipole-moment, polarizabilities nor chirality. The set of 
slow variables we consider are the standard slow hydro- 
dynamic variables of liquids, density n(q), current den- 
sity j(q), and temperature 6(q) (connected to energy 
conservation), and the slow structural relaxation of su- 
percooled liquids enters via a few memory functions. We 
thus provide a general framework for the analysis of light 
scattering spectra in supercooled molecular liquids which 
we expect will prove useful either for phenomenological 
discussions using fit-functions for the memory kernels — 
we will list all restrictions on these fit-functions — or 
for the consideration of specific scattering mechanisms. 
Our central technique for simplifying the spectra consists 
in small wave vector expansions of the memory func- 
tions as should be appropriate for disordered systems. 
Thus we adopt the idea of generalized hydrodynamics, 
which extends the regular hydrodynamic approach to 
larger frequencies. In detail, we use the one suggested 
by Gotze and Latz (T^| as it provides a physically reason- 
able description of glassy systems. Finally, we also derive 
Green-Kubo formulae which enable e. g. computer sim- 
ulations to determine the memory functions and thus the 
complete spectra directly. 

Our assumptions pertain to the systems under consid- 
eration and can be tested experimentally: We consider 
only the lowest non-trivial orders in wave vector transfer 
q, in order to find general results for the light scattering 
from amorphous dielectric and macroscopically isotropic 
and optically inactive materials within the framework of 
linear response, classical statistical mechanics and clas- 
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sical electromagnetism. The condition ga < 1, where a 
denotes either a typical molecular size, the average par- 
ticle distance or a collective correlation length, appears 
well satisfied for supercooled molecular liquids but ex- 
cludes studies of critical phenomena. Electromagnetic 
retardation also can be neglected for lu <C cq, where c is 
the speed of light. 

The general formulae for the spectra and con stitutive 
equations are presented in Section ||. Section III lists the 
central results which are discussed in Section IV , The re- 



sults for the depolarized spectra are compared to previous 
theoretical approaches in Section and conclusions in 
Section VI summarize our results. More technical aspects 



are contained in Appendizes A, B and C, and Appendix 
D outlines the application of our general results to spe- 
cific light scattering mechanisms. 



II. GENERAL FORMULAE 



A. Dielectric fluctuations 



The dielectric permeability fluctuates in space and 
time around its average and additionally a (diffusive) 
scattering occurs. In general, this scattering spectrum 
reflects the dynamical processes in the sample and de- 
pends on the frequency of the incident, u>i, as well as 
the scattered wave, u>/. A simplifcation is possible if 
one considers only small frequency shifts, u> — loi — u>f 
with \u>\ <C tOi- Then dynamical dielectric fluctua- 
tions (e^(q, i)|efei(q)) determine the scattering cross sec- 
tions completely M. Here the Kubo scalar product, 
(A(t)\B) = (l/k B T)(6A*(t)5B), is used, with T tem- 
perature and k B Boltzmann's constant. 

The fluctuation Seij (q, t) has even time reversal sym- 
metry, is a symmetric tensor of second rank, and, as 
Fourier-transform of a real quantity, fulfills q, t) = 
6tj(q, t)*. In particular, the long- wave-length limit 
£ij(q — > 0,t) is real. Different Cartesian components 
i,j, k, I are picked out depending on the polarization di- 
rections of the incoming and scattered light 0; see Ap- 
pendix A for more details. The dynamical evolution is 
given by the Liouvillian L via dtA = iCA. A Laplace 
transformation — convention f(z) = i dte tzt f(t) for 
3z > — thus leads to the problem to calculate, for 
q — > 0, the matrix elements of: 



(ey(q,z)|e fei (q)) = e„(q)| 
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kfcz(q) 



(i) 



The spectra at frequency to then are given by the 
imaginary part of (ey(q, w + i0)|e/-;(q)) denoted by 
(ey(q, cj)|efe;(q)) . They are real, positive and symmetric 
functions of frequency. 



B. Generalized hydrodynamics 

Light scattering measures large wavelength dielectric 
fluctuations. Even though qa <C 1 can thus be as- 
sumed, the limit q — > cannot be performed naively 
in Eq. ([!]). Because of density-, momentum-, and 
energy-conservation, there are poles in the resolvent 
R(z) = (C — z)~ x , which shift to vanishing frequency in 
this limit jlq] . Using the Zwanzig-Mori formalism, these 
hydrodynamic low frequency features can be identified. 
One introduces the reduced resolvent 



R'(z) = Q 
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QCQ 



-Q 



(2) 



where the projector Q = 1 — P projects perpendicular to 
the hydrodynamic modes: 



P = 



|n(q))(n(q)| , |0(q))(0(q)| 



(n(q)\n(q)) (9(q)\Q(q)) 



E 



|j«(q))(j-«(q)l 



(3) 

The standard hydrodynamic formulae are obtained in 
this approach if the reduced resolvent R'{z) is treated 
in a Markovian approximation, replacing its matrix ele- 
ments with frequency independent transport coefficients 
. Generalized hydrodynamics differs from this by re- 
taining the frequency dependence of R'(z) but still ne- 
glecting its wave vector dependence. This generalization 
is required for liquids at lower temperatures as the struc- 
tural relaxation slows down strongly. 

Following Ref. |Q we identify the fluctuating temper- 
ature 0(q) with the kinetic energy fluctuations, e K (q), 
that are orthogonal to the density fluctuations, Cy 0(q) = 
Q n e K {q) =e K {q)-n{<i)(e K \n)/(n\7i). Here c v = 3k B /2 
abbreviates the specific heat per particle of the kinetic 
degrees of freedom, and Q n is the projector orthogo- 
nal to the density. Conservation of the total energy 
e(q) = e K (q) + e p (q) implies 



c v CQ(q)=q^(q)-qj L (q] 



(gjn) 
(n\n) 



-£Qe p (q), (4) 



where superscripts L indicate the longitudinal part, j L = 
q - 5/q, eP (q) is the potential energy, and j e (q) the to- 
tal energy current. Note that since (0(q)|e p (q)) = 
one can replace Q n e p (q) = Qe p (q). The hydro- 
dynamic variables are orthogonal with normalizations: 
(n(q)|n(q)) = NS(q)/k B T, (j k (q)\ji(q)) = (N/m)S kl 
and (0(q)|0(q)) = NT/c v . Here, n is the average den- 
sity of N molecules, m the molecular mass, and S(q) 
is the equilibrium center-of-mass structure factor; Ski is 
Kronecker's symbol. Throughout the following we will 
neglect wave vector dependences caused by molecular 
length scales, and replace, e. g., the structure factor by 
its homogeneous limit given by the isothermal compress- 
ibility k t : S(q) = S(0) + O((qa) 2 ) with 5(0) = nk B Tn T . 

Considering, in generalized hydrodynamics ]lq] , the 
fluctuating temperature instead of energy fluctuations 
rests upon the experimental observation that the heat 
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conduction of glasses and liquids is not drastically dif- 
ferent. This aspect is discussed in Refs. where 
the generalized hydrodynamics also is tested by molec- 
ular dynamics simulations. This formulation of general- 
ized hydrodynamics accounts straightforwardly for a fre- 
quency dependent isochoric specific heat. 



C. Decomposition of dielectric fluctuations 

The exact resolvent calculus sketched in the previous 
section thus provides a reformulation of Eq. (0) , see Eqn. 
(A14) in Ref. Q. The reduced dynamics R'(z) and the 
projector P, which projects onto the hydrodynamic vari- 
ables, appear: 

(e«(q)|iJ(*)|e«(q)) - Mq)|P'(z)Mq)) + 
(ey(q)|(l - R'{z)C)PR(z)P{l - CR' (z))\e kl (q)) (5) 

Thus the coupling of the dielectric fluctuations to the hy- 
drodynamic variables is found; explicitly it is given when 
writing out PR(z)P in Eq. (||). Additionally there is a 
background spectrum, the first term on the right hand 
side of Eq. (§). 

Since the hydrodynamic modes have been projected 
out, (generalized) hydrodynamics postulates that the 
limit q — > in R'(z) can now be performed safely. This 
leads to the well-known result for Raman scattering [Q : 
The background spectrum consists of scalar, and sym- 
metric traceless tensor scattering 

(e ij (c i )\R'(z)\e kl (c l ))=S(z)6 ij S kl + 
T{z)(S tk 5 3l + S tl 6 jk - + 0(q 2 ) , (6) 

Explicit expression for S(z) and T(z) can be obtained 
by choosing special linear combinations of the dielectric 
tensor. Let soo = [e X x + £yy + £zz]/3 denote the long- 
wave-length limit of the scalar part and t 2 o — [2e Z z — 
£xx — £yy]/Vl2 the v = component of the correspond- 
ing spherical tensor t<x v \ (the prefactors are conventional). 
Then, S(z) = (s 00 \R' (z)\s 00 ) and T(z) = {t 20 \R' (z)\t 20 ) 
are the scalar and the tensorial correlations; see Ap- 
pendix B for a discussion of the applicability of these con- 
siderations to polar fluids. Both spectral contributions 
S(z), T(z) are autocorrelation functions of real variables 
with even time inversion symmetry. Their spectra thus 
are even and non-negative. 

In order to simplify the discussion of the couplings to 
the hydrodynamic variables in Eq. (0), it is useful to 
consider the corresponding generalized constitutive equa- 
tions [|l4],[l5|]. These describe the temporal decay of the 
deviations in a variable, say in <Jcjj(q, t), produced by 
an adiabatic perturbation with external fields coupling 
to the conserved variables, after the perturbations are 
switched of at time t = 0: 
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(&y(q,t)> =J2 ( SA «(t)) (A a \e tJ (q))/(A a \A a )- 

a 

/ dr(5A a {T)) {A a \iCR'{t~T)\e l3 {ci))/{A a \A a ) , (7) 
Jo 

where the (orthogonal) hydrodynamic variables are ab- 
breviated: Ax = »(q),A 2 = Q(q),A 3 = ^(q),^ = 
j y (q),A 5 = jz(q))- The special choice of the exter- 
nal perturbation considered in the constitutive equation, 
Eq. (0), prepares a fluctuation in (Seij(q,t)) which de- 
cays slowly because it requires decay of fluctuations of 
the conserved variables at large wavelengths, 0(l/q). 
Thus the time evolution of (5eij(q,t)) in Eq. (|?j), is 
determined by time-dependent couplings to the gener- 
alized hydrodynamics of the distinguished variables. As 
seen trivially via a Laplace-transformation, the identical 
couplings appear in Eq. (||) as in Eq. (Q), expressing 
that via these couplings, close to equilibrium, dielectric 
fluctuations at long wavelengths acquire slow hydrody- 
namic components. Equation (0) further identifies the 
non-hydrodynamic components contributing to the back- 
ground. The first term on the right hand side of Eq. 
(0) describes static or instantaneous couplings, whereas 
the second term describes dynamic couplings which need 
time to build up and may be characterized by a finite 
response time. These will become important when ap- 
proaching the glass transition upon cooling as then the 
structural relaxation times increase. 

Density fluctuations are coupled to the dielectric fluc- 
tuations statically via the scalar scattering mechanism 
only, as follows from Eq. (0) when inserting the projec- 
tor P from Eq. (§): 

(n(q)|(l - £#(*)) Mq)) = (n|s 00 )^ + 0{q 2 ) . (8) 

The dynamic coupling vanishes because Cn(q) = q k j k {<\) 
is again a hydrodynamic variable. Furthermore the scalar 
density cannot couple to the dielectric tensor fluctuations 
t 2v in the limit q — > 0. 

Similar arguments hold for the coupling of the tem- 
perature to the dielectric fluctuations. Since the kinetic 
energy is not conserved, there is a dynamic coupling in 
addition to the static one. In order to guarantee con- 
servation of the total energy Eq. (Q) is used. Observ- 
ing (£Qe p \R'(z) = (Qe p \ + z(e p \R'(z) and rearranging 
terms one finds to lowest order in q. 

c y (e( q )|(i-£i?'(z))| eiJ ( q )) = 

S tJ [(Q n e\s 00 ) + z(e p \R'(z)\ Soo )} + 0(q) . (9) 

Further couplings of order q can be ignored. The dy- 
namic coupling in Eq. (JsJ) arises from the separation of 
the total energy fluctuations into fast kinetic and (pos- 
sibly) slow potential ones. Thus the simplified handling 
of the reduced resolvent in this generalized hydrodynam- 
ics has to be paid by an additional frequency dependent 
coupling to temperature fluctuations. 

Finally for the coupling of dielectric fluctuations to cur- 
rent fluctuations 
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(j k (c i )\(l-CR'(z))\e ij (c l )) = 
~ El qi(T kl (q)\R'(z)\e tJ (q))/m + 0(q 2 ) , (10) 

one finds a purely dynamical coupling, as first recog- 
nized within so-called two variable models |2|, [| . A static 
coupling is excluded, since the currents jk and the di- 
electric tensor e;j have different time reversal symmetry. 
The dynamic coupling to the stress tensor, tm, which 
appears because of momentum conservation, £jfe(q) = 
Ei1i T ki( c l)/ m i can be evaluated in the long-wavelength 
limit: 

(r w (q)| J R'(z)|e ij (q)) = (p \R'(z)\soo)S ij 5 k i + 
{T2 \R'(z)\t 2 o)(5ikSji + 5 u 5 jk - h^u) + 0{q 2 ) ■ (11) 

Here p = [t xx + r yy + t zz }/3 and r 20 = [2t zz - t xx - 
T yy \/y/V2 denote the long- wavelength limits of the pres- 
sure and of the transversal stress tensor. 

According to the basic assumption of (generalized) hy- 
drodynamics jl4|,[ll|, one has to keep terms only to the 
order indicated, the remaining ones are assumed to be 
regular with respect to frequency z in the limit q — > 0. 

III. RESULTS 

From Eqs. (||) to ([ll]) the light scattering spectra fol- 
low if a scattering geometry is chosen and the appropriate 
tensor elements are calculated; see Appendix A for the 
used geometry. Polarizations vertical to, V, or in the 
scattering plane, H, are considered, where the standard 
abbreviation Ii a denotes polarizations for incoming and 
outgoing light. The spectra depend on q, lj and scat- 
tering angle 9. We find I^y(g,6',w) = IvHja,9,uj) as 
predicted by Rayleigh's reciprocity theorem [Q. 



the spectral shapes only, i.e. cause low lying hydrody- 
namic lines, but do not lead to long-ranged static correla- 
tions. Concurrently, the static couplings of the conserved 
variables to the dielectric fluctuations in Eqs. (|J) to (|Io| ) 
are g-independent for small wave numbers. See Section 
IV, why these results appear violated when considering 
the hydrodynamic limits. 



B. Depolarized spectrum 

Using the decomposition of the offdiagonal dielectric 
fluctuations, Eqs. (fjj), ( |To| ) and (|i"T|), one can identify 
a dynamic coefficient, ay#(z), which describes the cou- 
pling to the generalized hydrodynamic variables: 

a V H(t) = jj(T 20 \R'(t)\t M ). (15) 

It is a generalized elasto-optic or Pockels' constant fa- 
miliar from light scattering in solids 0, and is a real 
and symmetric function of time, as the two tensor vari- 
ables determining it are real with even time parity. Its 
Laplace transform, ayniz), therefore has an even spec- 
trum, a'y H (uj). This Pockels' constant describes the dy- 
namic coupling of the transverse current into the dielec- 
tric fluctuations and the constitutive equation becomes: 

9 

(5e V H(<l,t)) = -iqcos- J dr a VH (t - r)(<5j y (q, r)) . 

(16) 

From Eqs. (|||]) thus follows the general result for the 
depolarized spectrum: 



A. Total scattered intensities 

The total scattered intensities (except for standard 
prefactors |lj) can be obtained directly from Eq. (A4) 
in Appendix A and consist of scalar and tensorial scat- 
tering 0: 



Ivv(q,9) = (s 00 |soo) + g (*2o|t2o) 



lvH{q,9) = (t 2 o|*2o) 



lHH(q, 9) = cos 2 9 (soo|s o) + 
(l + \cos 2 9){t 20 \t 2Q ) . 



(12) 



(13) 



(14) 



The intensities are wave vector independent as the limit 
qa <C 1 is considered in systems where all molecular cor- 
relations are short-ranged. The conservation laws affect 



(ev/ff(q, z)\evH(q)) = T(z 

2 

q cos 



J 1 a VH {z) 2 CjAq,z) 



(17) 



Here Cj^z) = (j T (q, z)\j T (q)) 



denotes the correla- 
tion function of the transversal current fluctuations. The 
spectrum consists of a background arising from the sym- 
metric scattering in Eq. (|^), which commonly is dis- 
cussed as Raman line, and of couplings to the current 
fluctuations. Naively evaluating the depolarized spec- 
trum at vanishing wave vector would neglect this addi- 
tional contribution. It is small, of order 0(q 2 ), but is 
characterized by a time scale that diverges in the hydro- 
dynamic limit, q — > 0, and therefore dominates the low 
frequencyspectrum. The full current correlators appear 
in Eq. (|17|), which stresses that no assumptions about 
translational-rotational coupling are required in order to 
derive Eq. (|l7|). Explicitly this has been shown by the 
derivation of Eq. (|l7]) for a liquid of spherical particles 
in Ref . flSfl , which was tested in a simulation |20| . 
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C. Polarized Spectra 



For the VV spectrum, Eq. (|llf) suggests to introduce 
the elasto-optical constant 



avv(t) = o a vH(t) 



(p\R'(t)\s Q0 ) 
N 



(18) 



It has identical properties as ay nit) since the tensor vari- 
ables entering its definition again are real with even time 
parity. Another time dependent coupling function arises 
from the temperature fluctuations as described in Eq. 



ds m 
dT 



+ z 



(e p \R'(z)\s Q0 ) 
NT 



(19) 



where the thermodynamic derivative NT(dsoo/dT) n = 
(Q„e|soo) is written, which contains the total energy per- 
pendicular to density fluctuations fH| , Clearly, Eq. jl9| ) 
presents a generalized time or frequency dependent ther- 
modynamic derivative. The scalar variables entering its 
time or frequency dependent term are real with even time 
parity, and £,"(lv)/uj consequently is an even function of 
lu. The constitutive equation coupling the fluctuating 
hydrodynamic variables into 5eyv(q, t), from Eq. (Q), 
becomes: 



(5e V v{<l,z)) = (ds o/dn)T(Sn(q,z)) + 
e(z)(<59(q,z)) +a vv (z)q(Sj L (q,z)) , 



(20) 



where the thermodynamic relation (w| s oo) = 
NnKT(dsoo/dn)T is used. Collecting the terms in Eq. (||) 
one obtains when applying the mass-conservation law, 
which gives q 2 Cj J J (q,z) = qzC^(q,z) = z 2 C nn (q,z) + 
zNnKT^qC^^q.z) = zC en (q,z), 

(ew(q, *)|ew(q)) = S(z) + 4T(z)/3 
+zay V (z)Nukt + 2(dsoo/dn)Tavv(z)NnKT 
+ [(ds o/dn) T + za V v(z)] 2 C nn {q, z) 
H(z) 2 C e0 (q,z) 

+2[(ds 00 /dn) T + za vv (z)]£{z)C n@ {q, z) . (21) 

Here C nn (q,z) denotes the density-density correlation 
function, Cee(?, z) the temperature-temperature corre- 
lation function, etc; see Appendix C. Equation ( ^l| ) is 
our principal result for the VV spectrum. It extends the 
conventional hydrodynamic spectra to arbitrary frequen- 
cies. The small-wave-vector singularities are encoded 
in the generalized hydrodynamic correlation functions, 
C a p{q 1 z) 1 which are determined by the true resolvent, 
R{z), and can thus e. g. be obtained from simulations. 

From the information on the depolarized and the po- 
larized spectrum also the Ihh spectrum can be obtained 
even though it is not a simple linear combination. The 
fluctuating variable coupling to the distinguished vari- 
ables for the HH scattering in the geometry of Appendix 
A is given by 



(Se HH (q, z)) = - cos 9 (ds 00 /dn) T (Sn(q, z)) 
-cos9 £{z){5B(q,z)) 

-[a HH {z) cos9 - a VH (z)}q(Sj L (q, z)) , (22) 

where we abbreviated clhh(z) = avv(z) — clvh(z)- The 
general spectrum in the HH geometry now reads 

{e HH (q, z)\e HH {qj) = S(z) cos 2 9 + T(z)(l + J cos 2 9) 
+z [o,hh(z) cos9 — avH(z)] 2 Nukt 

' cos 9 [auH (z) cos 9 — avn{z)] Nnux 



T 

ds 00 \ 



—j +za HH (z) 



H(z) 2 cos 2 9 C &@ (q,z) 



cos 9 - za VH (z) \ C nn (q,z) 



-2£(z) cos 6>{ 



ds 00 \ 
dn J , 



zclhh{z) 



cos 9 



-za V H 



{z)\c n e(q,z) 



(23) 



One notices that even for a general molecular fluid, trans- 
verse current fluctuations do not couple into the HH 
spectrum. 



D. Generalized Green— Kubo relations 

Ten frequency-dependent matrix elements built with 
the reduced resolvent R'(z) have been identified in the 
expressions for the light-scattering spectra in supercooled 
liquids. Five generalized transport coefficients and ther- 
modynamic derivatives are needed in order to describe 
the correlators of the hydrodynamic variables G3] . They 
are the shear viscosity K s (z), the thermal conductivity 
A(z), the dynamic specific heat cy(z), the dynamic ten- 
sion coefficient /3(z), and the longitudinal stress relax- 
ation kernel Ki(z) (explicit expression are summarized 
in Appendix C) . The remaining five frequency dependent 
kernels encode the details of the light-scattering process: 
S{z), T(z), clvh(z), avv(z), 

These expressions involving reduced resolvents are very 
suitable for approximations, since they do not exhibit 
hydrodynamic singularities. To determine them from 
other theoretical approaches or from computer simula- 
tions, however, it is more convenient to find a formula- 
tion in terms of correlation functions involving the full 
dynamics. For the considered cases at q = 0, this is 
made possible by the conservation laws which allow one 
to derive Green-Kubo relations for the memory functions 
(transport coefficients) expressing them in terms of (time 
integrals of) autocorrelation functions of the correspond- 
ing fluxes Jl5|,pL6||. From the identity @: 



(X\R'(z)\Y) = (X\R(z)\Y) + 

(X\(l - R / {z)C)PR(z)P(CR'(z) - 1)\Y) 



(24) 
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one observes that the reduced matrix elements can be 
rewritten as full matrix elements and correlation func- 
tions of the hydrodynamic variables contained in PR(z)P 
with frequency dependent coefficients. Since for q — > 
the coefficients involving Ln q and Lj q vanish due to par- 
ticle and momentum conservation, only the temperature 
fluctuations contribute to the frequency dependence of 
the coefficients. To derive generalized Green Kubo re- 
lation, we only need variables X — QX and Y = QY, 
respectively. Therefore all static couplings X\PR(z)P to 
the hydrodynamic variables in ( p4| ) will vanish too and 
the generalized Green-Kubo relation is given by 

(X\R'(z)\Y) = (QX\R(z)\QY) + 

[(X\Qe p ) + z(X\R\z)\c p )\[{Qc p \Y) + z(e p \R'(z)\Y)\ 



NTzcv(z) 



for q — > . 



(25) 



Here we made use of Eq. (|4|) and (C3) in the limit q — » 
and of the identity R'(z)LQ = Q + zR'(z). Rotational 
invariance implies that the second term in Eq. J25| ) is 
nonzero only for scalar variables X, Y. Thus, e.g. the 
standard Green-Kubo relations for the shear viscosity, 
X = Y = T20, and heat diffusion, X = Y = j e , are 
found. Moreover, it follows that in the elasto-optic con- 
stant clvh{z), Eq. ( [15]) , and in the tensor background 
spectrum T(z), Eq. (ro|), the reduced resolvent can be 
replaced with the full dynamics. 

In order to obtain tractable expressions for the 
remaining kernels and make contact with the stan- 
dard Kadanoff-Martin approach [fl8|| , we define an- 
other, the conventional, fluctuating temperature by 
cyT(q) = e(q) — n(q)(e|n)/(n|n), with normalization 
(T(q)|T(q)) = NT/cy. Also, let Q denote the projec- 
tor orthogonal to density, currents and T(q). Choosing 
X = Y = e p in Eq. @ one finds 



{Qe p {z)\Qe p ) c v (z)-c v (c v (z) - c° v ) 2 



NT 



zcy(z) 



(26) 



Since Qe p — Qe p + (cy — Cy)T the left-hand side im- 
plicitly contains hydrodynamic poles due to energy con- 
servation. However, 



(Qe p (z)\Qe p ) (Qe p (z)\Qe p ) (c v 



c° v ) 2 



NT 



NT 



ZCy 



(27) 



and, according to Kadanoff-Martin, the first term on the 
right hand side is free of poles in the limit z — > [ p"8[ . 
Combining Eqs. (|||27|) one derives 



c v (z) = c v 



„0 \2 



(c°y) 2 



zc v {Qe p (z)\Qe p )/(NT) 



(28) 



In the liquid phase the dynamic specific heat attains its 
thermodynamic value cy(z) — > cy for z — ► 0. Equation 
( p8| ) demonstrates explicitly that the Gotze-Latz resol- 
vent R'{z) indeed devoids all hydrodynamic singularities 



and is compatible with the conventional Kadanoff-Martin 
formalism. It differs in an explicit frequency dependence 
of cy(z), which arises from the splitting of the conven- 
tional temperature fluctuations, T(g,t) into fast, the ki- 
netic ones Q(q,t), and structural slow ones. 

Similarly, substituting X = p, Y = e p in Eq. (|25| ) 
yields 

(Qp(z)\Qe p ) _ [3(z)~[3 [/?(z)-/3°][c v (z)- Cv ] 



NmT 



zc v {z) 



, (29) 



where 0° = {p\Q n e K )/(NmT). Since Qp = Qp + 
m(3°Qe p /c v + m{f3 - (3°)f the left hand side can be 
written as 



NmT 

(3° 



(Qp(z)\Qe p ) 



1 



NmT 



(Qp(z)\Qe p ) 



NTc° v 



ZCy 



Again the memory kernels appearing on the right-hand 
side are regular in the low-frequency limit. Collecting 
terms leads to the generalized Green-Kubo formula for 
the dynamic tension coefficient 



Cy 



NmT 



(31) 



In a similar fashion the corresponding Green-Kubo for- 
mulae for the dynamic temperature coupling, the scalar 
background spectrum, the remaining term contributing 
to ayy(z), and the longitudinal stress-stress correlation 
function can be obtained: 

_ dsoo \ c v (z) c v (z) (Qs 0Q (z)\Qe p ) 



dT 



Cy 



S(z) 



NT£(z) 2 NT <9s 00 x 2 



ZCy(z) ZCy 8T ) n 

(p\R'(z)\s 0Q ) mf3(z)T 



NT 

— + (Qs 00 (z)\Qs 00 ) , (33) 

m(3T OsqqX 



N 



+ 



ZCy(z) ZCy Ol 

(Qp(z)\Qs 00 ) 



N 



Ki{z) 



mT , , mT 

— p{zf /3 2 

ZCy(Z) ZCy 



(Qr L (z)\Qr L ) 
Nm 



(34) 



(35) 



In particular, the last term on the right-hand side of 
Eq. ( |35| ) is related to the the longitudinal viscosity 
ir)l = (Qt L ( z -* 0)\QT L )n/N. Therefore in the low- 
frequency limit where cy(z) = cy +izc' v , (3{z) = f3 + iz(3" 
one finds 



Ki{z^Q) 



2imTf3 



0" 



cy 



(36) 



To summarize, all ten frequency-dependent kernels can 
be expressed in terms of the full resolvent R(z), and 
therefore they can directly be obtained from molecular- 
dynamics simulations. 



6 



IV. DISCUSSION 



A. Depolarized spectra 



The depolarized frequency dependent spectrum, Eq. 
(|l7|), shall be discussed in detail as it provides the most 
compact expression but also exhibits clear qualitative 
changes when supercooling the liquid. It consists of three 
frequency dependent contributions, a background, the 
Pockels' constant and the transverse current correlator. 

The current correlators can be taken from theories for 
the dynamics of the liquid under study or from com- 
puter simulations. Alternatively, the generalized hydro- 
dynamic approach shifts the problem to calculate the 
transverse correlator to the problem to calculate correla- 
tions of the transversal stress tensor, viz. the frequency 
dependent shear modulus, K s {z) — {t 2 o\R'(z)\t2o) n/N, 
where the resolvent devoid of hydrodynamic fluctuations 
from Eq. (J2J) appears again. The separation of the hydro- 
dynamic poles from structural relaxation thus is achieved 
leading to Q : 



CjAq,z) 



-(N/m) 



z + q 2 K s (z)/nm 



(37) 



The result from hydrodynamic theory for the depolar- 
ized light scattering from equilibrium molecular liquids 
can be obtained if the frequency dependence of mem- 
ory functions built with the reduced resolvent R'{z) is 
neglected by using the Markovian low frequency limit. 
Then the depolarized Pockels' constant becomes purely 
imaginary: 



sra(z) — > ia' 



VH 



dt a VH (t) for z -> . (38) 



Therefore, and using the standard hydrodynamic result 
for the shear viscosity, K s (z — > 0) — > ir\ s , one recovers the 
result for the depolarized spectrum first obtained within 
simple models in Refs. [ 



I VH {q,e,u) hy . h =T"{u> = 0) 

2 ,» )2 2 Nq 2 Vs /(nm 2 ) 

\ a vH) co& 7; 2 1 1 1 — / — \2 

2 uj z + (q z rj s /nm) z 



(39) 



The transverse momentum diffusion cuts a central line 
of half-width q 2 r] s /(mn) and amplitude proportional to 
{a'y H cos ^) 2 nN/rj s out of a flat background. Note, that 
within hydrodynamics the background is structureless, 
and that the spectrum is positive owing to an elemen- 
tary Schwartz inequality, (^ol^oX^ol^o) > (t 2 o|^2o) 2 - 

The structural relaxation of liquids cooled down to and 
below the melting temperature slows down strongly and 
the frequency dependence of the memory functions with 
reduced resolvent cannot be neglected anymore; for a dis- 
cussion of structural relaxation see e. g. the review 
The result in Eq. (|l7]) can handle this situation, as the 
memory functions built with R'(z) — there are 3 in Eqs. 



(|T^,^) — may either be modeled appropriately or can be 
taken from other theories or simulations. Within the gen- 
eralized hydrodynamic approach a glass, or amorphous 
solid, is obtained whenever a structural relaxation pro- 
cess, with time scale r Q , is slow compared to the hydrody- 
namic frequencies. Assuming further that the dynamics 
in R'{t) at shorter times, denoted by 773, can be neglected 
for this frequency range, then the memory functions in 
Eq. ( |l7j ) can be approximated by: 



K s (z) 



iT. q 



a VH (z) 



"VH 



(40) 



for l/r Q <C \z\ <C I/773. Note, that this is equivalent 
to time-independent values, K s (t) — Goo and flyg(t) = 
a,y H for r^«i <C T a . Therefore, the poles in Eq. (HO) 
are called non-ergodicity poles as they describe frozen- 
in, non-relaxing components. = mnc 2 -, is the glassy 
shear modulus familiar from Maxwell's model and a^ v 
is the Pockels' constant (often denoted V44) quantifying 
the elasto-optic coupling in the glass 0. Whereas G^o 
and r s need to be positive, the sign of a^ v is undeter- 
mined; an imaginary part in anv exists in principle but 
does not contibute to the spectrum in the hydrodynamic 
limit. Equations ([l7],[l0|) predict for the hydrodynamic 
glass spectrum: 



I VH (q,6, u)hy.g, = Tg + 

fa cos -a 00 ) 2 q 2 Ts/(mn) 

(q cos 2 a VH ) ^ 2 _ q2( a )2 + {ujq 2 Ts/nm) 2 



(41) 



Two transverse phonon peaks characterized by the trans- 
verse sound velocity, Ct, and a width cx g 2 r s , appear 
which are described as damped harmonic oscillations. 
The background consists of a central line which can- 
not be resolved and a structureless continuum, T{z) — 
-T^/z + iTg". 

Note, that both hydrodynamic expressions, Eqs. 
(39,fy]), violate the sum-rule for the total intensity, Eq. 
(13), and predict wave vector dependent total scattered 
intensities. The reason, of course, are the Markovian ap- 
proximations in Eqs. ( |l7| , [40"| ) which fail to describe the 
dynamics outside the hydrodynamic range. Non-trivial 
spectra obtained in glasses on frequency scales charac- 
terized by 773 p2[ also require more elaborate expressions 
for the memory functions in Eq. (fl7|). 

For temperatures around the liquid-to-glass crossover 
at T c , neither the assumption r a uj 3> 1, nor the estimate 
Tpuj <C 1 hold and the depolarized spectra exhibit anoma- 
lies |2|| . The mode coupling theory of the structural re- 
laxation there suggests to model the reduced resolvent 
as K s (z) « -Goo[l - (1 - izTa)- fiaD + (-izt Q ) a ]/z for 
T w T c , and similar expressions for the other two mem- 
ory functions. Whereas the Cole-Davidson behavior (the 
first two terms) is a (rough) model of the a-process, and 
Pcd will differ for different resolvent matrix elements, the 
power law with exponent a describes the universal "criti- 
cal" decay close to the transition. Here to is a microscopic 
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time and the exponent a and the (true, universal) expo- 
nent b of the high frequency von Schweidler wing of the 
a-process, K s (l/r a <C z <C 1/tb) ~ (—izT a )~ b /z, are 
related; see e. g. the review for further information. 



B. Polarized spectra 

A most prominent feature of the polarized spectrum is 
the Brillouin peak arising from propagating sound waves. 
Upon cooling the liquid, structural relaxation manifests 
itself predominately by a gradual change of the sound ve- 
locity and the damping constant. Considering the enor- 
mous increase of the transport coefficients, e. g. the lon- 
gitudinal viscosity which describes the damping in the 
liquid, this clearly points out the necessity to consider the 
frequency dependence of the reduced resolvent as done in 
generalized hydrodynamcis. Furthermore, the couplings 
avv(z), £(z) as well as the background spectra S(z), T{z) 
exhibit non-trivial z-dependences in the frequency regime 
of interest. 

The structural relaxation of the Pockels' constant 
ayv{z) (with numerically different constants) can be 
modelled as given in Eqs. (|3|) and ( f4(i| ) and as described 
at the end of the previous section. The explicit factor 
z in the frequency dependent part in £(z), Eq. (19), 
cancels a possible non-ergodicity pole. The coupling to 
temperature fluctuations interpolates smoothly between 
its low-frequency thermodynamic value (dsoo/dT) n and 
a high-frequency coupling, £00, characteristic for a glass. 
A renormalization also appears in the effective coupling 
to the density fluctuations, which in Eq. ( pl|) is described 
by the Pockels' constant V12 = (dsoo/dn)T + zavv(z) 
and in Eq. ( |23| ) by (V%2 cos 6 — V44), respectively, where 
P44 — zcivh(z). Note, that the frequency dependent 
renormalization of (dsoo/dT) n in Eq. ( |l9| ) vanishes, if 
only the hydrodynamic fluctuations, density and tem- 
perature, contribute to the scalar scattering; then also 
avv(z) = 1a HV (z). 

In order to obtain the spectrum in the true hydrody- 
namic limit one substitutes the appropriate correlation 
functions — see Appendix C — and replaces all memory 
functions with their low frequency limits; only in cv{z), 
(3{z) and £(z) linear terms in z need to be kept as can be 
seen from Eqs. (|2^|3l],[32|) . There are the three familiar 
hydrodynamic resonances superimposed on the Raman 
background: the Brillouin doublet of sound modes and 
the Rayleigh heat pole. The spectrum is obtained from 
determining the residues of these poles to lowest order in 
frequency and wave vector. 

I vv (q, 9, io) hy . L = (Q Sqo (lo = 0)|Qs 0O )" + 4T"(w = 0)/3 + 

rX Br - cVr, 



NuktVi 



+x 



R.7- 



1 



7 (w 2 - 
q 2 D 7 



c 2 q 2 ) 2 



M 2 r ; ) 2 



where the adiabatic sound velocity c = y/ 7 / JrnnKr) , the 
longitudinal sound damping Ti = £>t(7 — 1) + ry//(mn), 
and the heat diffusion constant Dp = A" ' /cp appear (see 
Appendix C for details). The Pockels constant is given 
by the thermodynamic derivative, V12 — dsQo/dn)T, and 
the flat background consists of scalar and tensor parts. 
Neglecting contributions from temperature fluctuations, 
X K - — X Br - = 1, one regains the well-known result 
for light scattering from (hydrodynamic) density fluctua- 
tions Then, the Landau-Placzek result (7 — 1) is 
recovered for the relative intensity of the Rayleigh to the 
Brillouin lines, where 7 = cp/cy denotes the ratio of the 
isobaric to isochoric heat capacity p], [?|, |l5|] . In the gen- 
eral case, scattering from temperature fluctuations leads 

2 

to (presumably small) corrections: X Ri = (1 — -f-p^) 2 

and X Bt ' = (1 + (7 - 1)|^) 2 , where £ = ds m /dT) n . 

For the glass, ut q 3> 1, actually the identical formula 
holds where, however, the isothermal compressibility, the 
sound velocity and damping constants, as well as the 
couplings and the background — compare Eqs. (39 4l]) 
- are renormalizcd. The high frequency values of the 
memory functions appear in the formally identical def- 
initions of c, Dp and T;, where, in the frequency win- 
dow 1/t q <C lo <^ I/t/3, simple Markovian expressions 



7 L0 2 + (q 2 D T ) 2 



(42) 



like cv{z) 



izc 



V.g 



are appropriate for ergodic ma- 



trix elements (cy(z), /3(z), £(z) and A(z)), and frozen-in 
components, leading to Ki(z) — —K^/z + iK'/ g , ap- 
pear in Ki(z), ayv(z) and in the Raman background 
lines. Thus, e. g. the expression for the Pockels con- 



stant becomes: = ds o/dn)T 



l vv- 



A non-trivial 



renormalization appears in the isothermal sound velocity 
or cquivalcntly the isothermal compressibility, because of 
the frozen structural relaxation in the longitudinal fric- 
tion function: 



oo\2 



1 



ran 



(43) 



as first observed by Mountain || and predicted from mi- 
croscopic expressions by the mode coupling theory pT| . 
Note that the reduction of the Brillouin and Rayleigh 
intensities described by Eqs. (^,^3|) is caused by the 
appearance of frozen-in density fluctuations (an elas- 
tic Mountain-line) which contributes the missing weight, 
NnKT((ds m /dnf T -(j§,f{V^) 2 ). 

T 

For the HH spectrum let us just mention that for 
9 = n/2, Eq. (2^) simplifes considerably, since the 
scalar scattering drops out completely. For frequencies 
uT a <^ 1, avn{z) —> ia!y H holds again, and one finds 

I H H{q,e= ^,Lo) hy . L = T"(u = {))- (uja' VH ) 2 C nn {q^)" 

(44) 

The hydrodynamic modes are suppressed by a factor of 
1 ■ the Brillouin line cuts out a Lorentzian with 



e.g. 



half-width q Ti and amplitude N (a! v H ) / (mY 1) out of a 
flat background. 
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C. Intensity ratios 

Interesting intensity ratios can be constructed if the 
scattered intensities with different polarizations of the 
light before and after the scattering process are consid- 
ered. The standard depolarization ratio compares the 
off-resonance (u> 3> cq) frequency-dependent background 
intensities neglecting all hydrodynamic modes: 

IvH(q,6,v) = T>) < 

l vv {q,0,u) (Q S0 o(w)|Qs o)" + 4T"(w)/3 _ 

(45) 

The expected depolarization ratio 3/4 is recovered if the 
only scalar scattering mechanisms are density and tem- 
perature fluctuations. 

V. COMPARISON WITH 
PHENOMENOLOGICAL APPROACHES 

It appears worthwile to discuss earlier phenomenolog- 
ical approaches within our general framework. Again wc 
shall concentrate on the depolarized spectrum which has 
been the focus of a long list of theoretical descriptions. 

The original approach of Andersen and Pecora ||] 
and Keyes and Kivelson || captures the spectra as pre- 
dicted by hydrodynamic theory. Earlier approaches had 
failed to do so, because they overly constraint the hy- 
drodynamic transport coefficients. The two or three 
variable approaches in go beyond hydrodynamics 

as they model the background spectrum as a (sum of) 
Lorentzian(s). But they miss the frequency dependence 
of the shear modulus and of the Pockels constants and 
thus do not present a valid extension beyond the hydro- 
dynamic result in Eq. (p9[). Consequently their results 
could not be applied to supercooled liquids. They also 
failed to recover the sum rule in Eq. (|13|). Vaucamps et 
al. extended this approach by using Maxwell's model, i. 
e. the approximation K s (z) = —mncj,/[z + i/T s ] for the 
shear modulus p|. More elaborate models, which iden- 
tify slower and faster relaxational processes in the shear 
modulus K s (t), were sugg ested in Refs. [@-[Ll). Models 
with more than two variables can be brought into a form 
suggested by a simple viscoelastic approximation (VA) 
to our general result, Eq. (|l7|): 

$%(z)<x ~* + 7 - = rr-Ti ( 46 ) 

VHK ' z + iTr [z + ij VH ] 2 z-q 2 c^/[z + i/T s ] y ' 

where A = {qa^ H cos f) 2 JV/(mT°°) and Eq. © 
gives a simplified a-process only description akin to 
Maxwell's model. The multi-variable phenomenologi- 
cal models, however, rest on specific assumptions about 
the rotational-translational coupling, and are, therfore, 
forced to introduce constraints on the coupling parame- 
ters in order to recover Eq. (fig); see Ref. |R| for a de- 
tailed discussion of this aspect. Note that the a-process 



only model, Eq. ([46]) , fails to describe higher frequencies 
as can be seen from its incorrect prediction of a vanishing 
damping of the shear waves in the glass, see Eq. fl4l|), 
and from its failure to recover the exact sum rule, Eq. 

©■ 

The introduction of retardation effects via memory 
functions in the work of Wang [jlO| and Dreyfus et al. 

removed constraints on the coupling parameters, 
and, within the latter phenomenological approach, an ex- 
pression formally equivalent to Eq. (Il7]) was given. 



VI. CONCLUSION 

In this paper we discussed the light scattering spectra 
for a one-component molecular liquid incorporating slow 
structural relaxation and thus extending the description 
to supercooled liquids and glasses. In contrast to earlier 
phenomenological approaches no assumptions are made 
on the nature of the scattering mechanisms and on the 
origin of the structural relaxation. Our results are akin to 
the theory of neutron scattering experiments where it is 
shown that, in general, density fluctuation functions are 
measured without the need for special models of their 
dynamics |l^,^6). For example, Eqs. ( p7|j3^ ) indicate 
how the shear viscosity can be measured by light scat- 
tering, without specifying what microscopic dynamical 
mechanisms contribute to the decay of transversal cur- 
rents. In particular, we clarify the question of how many 
generalized frequency dependent elasto-optical constants 
and memory functions are needed for a general descrip- 
tion of the long-wave-length fluctuations of the dielec- 
tric tensor. Furthermore we give explicit microscopic 
formulae for the background spectrum and the elasto- 
optical constants which could serve as a starting point for 
approximations once a choice for the scattering mecha- 
nism is made. For isotropic particles considering dipolc 
induced-dipole scattering this has been performed in p9| . 
Alternatively, computer simulations could be employed. 
Our central assumption is to consider systems charac- 
terized by short-ranged equilibrium correlations. Thus 
we can neglect the wave vector dependence in matrix 
elements describing the structural relaxation leading to 
Green-Kubo like formulae for the memory functions as 
familiar from simple hydrodynamics. 

Whereas previous approaches were mainly concerned 
with the nature of the depolarized spectrum and corre- 
spondingly e. g. attribute the scalar fluctuations to a 
combination of density and temperature fluctuations, we 
carefully make the distinction of dielectric fluctuations 
coupled to hydrodynamic modes and the ones orthogo- 
nal thereto. Consequently one obtains a theory that com- 
bines conventional Brillouin and Raman scattering. For 
example, molecular vibrations or rotational motion give 
rise to scalar as well as tensor scattering which appears 
as background to the hydrodynamic resonances and is 
included in our framework. 
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An important aspect of our results is that no assump- 
tions on e. g. translational rotational coupling or about 
the concrete description of the structural relaxation were 
necessary. The general aspects of the light scattering 
spectra were worked out and correlation functions were 
defined, which can, in principle at least, be measured ex- 
perimentally, and which contain the general information 
about the dynamics of the sample under study. 
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APPENDIX A: SCATTERING GEOMETRIES 

The scattering plane, i.e., the plane that contains both 
the wave vectors , ky of the incident and scattered 
wave, is chosen as the x—z plane. The direction of the 
momentum transfer q = ki — k/ is taken to be (0, 0, —q). 
Since we consider only small frequency shifts |kj| ~ |k/| 
the scattering angle 9 — Z(k/,k<) is related to the mo- 
mentum transfer via q = 2ki sin |. The dielectric fluctu- 
ations corresponding to the conventional scattering ge- 
ometries, viz. polarization perpendicular (V) to and in 
the scattering plane (H) , are then given by 



Se vv (q) = Se yy (q) 



5evn(q) = 5e xy (q) sin - - 5e yz (q) cos 

,., flH j iu - -r u€ y 



2 "f-v-v — 2 
9 

5e H v(q) = Se xy (q) sin - + 6e yz (q) cos • 



<5e/m(q) = Se xx (q) sin 2 - - 5e zz (q) cos 2 - (Al) 

The scattering intensities are expressed via the spectra 
of the dielectric fluctuations 

I io (q,6,uj) = (e io (q,w)|e io (q))" 

i 

dtcos(ujt)(e io (q,t)\e io (q)) (A2) 

The scattering intensities obey the sum rule 

±- f ckjl io (q,6,w)=l io {q,6) (A3) 
27r Jo 

where the total intensities ho(q, 9) are determined by the 
thermal fluctuation of the dielectric tensors 

I lo (q, 9) = (e io (q, t = 0)|e io (q)) = ^f^l . (A4) 



APPENDIX B: POLAR FLUIDS 

Due to the long-ranged character of dipolar interac- 
tions, fluids of molecules possessing a permanent (elec- 
tric or magnetic) dipole moment need additional consid- 
erations. Surface polarizations (in the case of electric 
dipole moments) or demagnetizing fields (for magnetic 
moments, where we will use the electric language in the 
following) arise for an uniform sample and lead to a di- 
electric response to the applied external (vacuum) field 
which depends on the shape of the sample |l6|,|24|]. Thus 
our use of rotational isotropy in the main text appears 
not justified. However, because of the strong angular de- 
pendence of the dipolar interaction, polar fluids can lower 
their free energy by the formation of domains or textures, 
and macroscopically isotropic response functions and di- 
electric fluctuations can be expected [P^-p6|. For these 
situations, where ferroelectric order is absent, and where 
"chain-formation" |^7j does not lead to large aggregates, 
i. e. where qa <C 1 holds, our general formulae apply. 



APPENDIX C: GENERALIZED 
HYDRODYNAMICS 

The correlation functions in t he g eneralized hydrody- 
namic limit can be expressed as [ |14[ 

C nn (q,z) = ~Nnn T /[z- (c q) 2 /[z + q 2 K t (z) / (ran) 

~q 2 mTf3(z) 2 /[zc v (z) + q 2 X(z)]]] , (CI) 

C n @(q.z) = -NTq 2 /3(z)/[[z 2 - (c q) 2 

+zq 2 Ki{z)/{mn)][zc v {z) + q 2 \{z)\ - q 2 mTz(3{z) 2 ] , (C2) 

Cee(q, z) = -NT/[zc v (z) + q 2 X(z) 
-q 2 mTz(3(z) 2 /[z 2 - (c g) 2 + zq 2 K^z) / {mn)]\ . (C3) 

Here cq = (mnKT)" 1 ^ 2 denotes the isothermal sound ve- 
locity, Ki(z) = (T L \R'(z)\T L )n/N the longitudinal stress 
relaxation kernel, and X(z) = (j£ \R'(z)\j^)/ (NT) the 
thermal conductivity. The frequency-dependent tension 
coefficient and specific heat are given by 



f3(z) =(3 + z(p\R'(z)\e p )/(NmT) , 



c v (z) =c v + z(e r \R'(z)\e^)/(NT) 



(C4) 



(C5) 



where (5 = (p\Q n e)/(NmT), c v = (e\Q n e)/(NT) retain 
their standard values. 

There appear three hydrodynamic modes in the liquid 
state. First, the sound doublet z = ±cq — iq 2 Ti/2 + 0(q :i ) 
with the adiabatic sound velocity, c 2 = c§ + mTf} 2 /cy 
and 
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K 1 / 

nm 



mT(3 2 , \" 



C 2 Cy 



Cv 



(C6) 



where we have written iK" — Ki{z — * 0),iX" = X(z — > 
0), c v (z -> 0) = c v + izc' v , /3(z -> 0) = /3 + iz/3". In the 
generalized hydrodynamic approach the sound damping 
consists of four parts: viscous friction, losses at conver- 
sion from mechanical into thermal energy, losses from 
storing and extracting thermal energy, and thermal dif- 
fusion. Although the expression looks unfamiliar, the as- 
sociated Green-Kubo formula, Eq. (|36|), shows that the 
longitudinal viscosity appears, leading to the familiar re- 
sult 0,^5)- Second, there is the heat mode z = —iq 2 Dx, 
where the thermal diffusion constant reads 



D T = = 

C z Cy Cp 



(C7) 



For the second equality in (C7) we made use of 
standard thermodynamic transformation formulas. The 
merit of the generalized hydrodynamic approach lies in 
the separation of the long-wavelength and low-frequency 
properties |jl4| . The only wave- vector dependence arises 
from the conservation laws. The structural relaxation 
is captured in generalized transport and thermody- 
namic derivatives, which exhibit a significant frequency- 
dependence in the regime of interest. 



APPENDIX D: SPECIAL SCATTERING 
MECHANISMS 

In order to show that all contributions to the light 
scattering spectra discussed in the text will be present in 
general, and also in order to exemplify how our general 
formulae can be used, specific scattering mechanisms are 
listed: 

1. If single scattering processes from biaxial (likely chi- 
ral) molecules are considered, the background spectrum 
arises from the two irreducible spherical tensors: 



soo = -^(cixx + oiyy + a zz )n(q) , 



0, 



(Dl) 



2v = ±={a xx - ayy) £[I>$(nO + D^a(fii)] 

* i 

+ 77=(2a^ - °<- xx ~ a yy ) ^2 D^li^i) > 



where the scalar dynamics follows the density fluctua- 
tions and and the Wigner functions Dup (fi) capture the 
dynamics of the molecular orientation in terms of Euler 

angles (l = (<p t 4, X ) & 

2. In the special case of linear molecules, a xx = a yy , 
the tensorial fluctuations simplify 



t2v 



-(a zz - a xx ) ^2 Y 2v(0i, Wi) 



(D2) 



and the standard spherical harmonics appear M. 

3. In the case of optically isotropic particles, depolar- 
ized scattering can arise from first order dipole induced 
dipole (DID) scattering Q: 



tlv — 



d 3 k 



V2J (2 



»(-k)T 2 „(k)n(k) 



(D3) 



-4-7r-\/87r/15Y2i/(k) is the static dipole 



where T 2j ,(k) 
tensor. 

4. If second order DID scattering of isotropic particles 
is considered, then the scalar scattering consists of a con- 
tribution coupling to the density fluctuations, Eq. (Dl), 
and of another non-hydrodynamic contribution: 



a 

T 



E 



(2tt)6 ( } 



T 2 „(k)7Wp)n(-k- p)n(k)n(p) 



(D4) 



/i v — i-i —v 



r 2 , tf _ p (k)T 2iM (^)n(-k - p)n(k)n(p) . (D5) 



This implies that in this case the depolarization ratio 
from Eq. (B) is smaller than the expected value 3/4. 
Furthermore the scalar fluctuations overlap with the or- 
thogonalized energy, i.e. (soo|Qn e ) is nonzero, and the 
coupling of dielectric fluctuations to temperature are rel- 
evant. In particular, one finds that even in the hydrody- 
namic limit for a liquid the Landau-Placzek ration is not 
fulfilled. 

5. The lower depolarization ratio also already arises 
from first order DID of linear molecules: 



«oo = -^a(a zz - a xx ) x 

£(-!)" J 7^)3 T 2 ^(k)n(-k)n 2 ,(k) (D6) 



2 2 2 
fi v — fx —v 



(2tt) 3 

T 2 ,^(k)n(-k)n 2/t (k) , (D7) 



where the fluctuating tensor density reads n 2/J (k) = 

\f&Kf\5>Yji Y 2i*(0i,<Pi) exp(-ik -r t ). 

6. If the molecules are not considered as rigid there 
may be Raman active modes due to intramolecular vibra- 
tions. If Q a> i is the vibrational coordinate and (da/dQ a ) 
the corresponding change of the polarizability for mode a 
there is a contribution to the scalar dielectric fluctuations 
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s °° = EEao-^ (D8) 

i a 

Here the sums run over all molecules and Raman modes. 
A similar expression can be derived for the dielectric ten- 
sor fluctuations. 
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